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1. INTRODUCTION 

Since their discovery (Klebesadel, Strong, & Olson 1973) 
cosmic gamma-ray bursts (GRBs) have been the focus of 
observations made with numerous experiments. (For an 
extensive review of burst observations see Fishman & 
Meegan 1995). Due in large part to the lack of a distance 
sca!e the origin of GRBs is still unknown. Developing 
detailed burst models is quite challenging when the intrinsic 
energy budget in gamma rays is uncertain by ~ 20 orders of 
magnitude. Burst models must therefore rely on the obser- 

straints ^ * Vane * y ° f burSt P ro P erties for meaningful con- 

Among the observed burst properties that can prove to 
be useful are the timescales over which these transient 
emi l radiatlon - Man y burst models (see Dar et al 
qq? ; Sf* & Me f aros 1992 ; Usov 1992; Begelman et al. 
993; Meszaros & Rees 1993; Woosley 1993; and Katz 
1994) make predictions that should be observable in a GRB 
time history. The burst time histories may be shaped on a 
number of timescales by either the emission process (eg 
through cooling timescales), the burst environment 
(through interaction with intervening matter, beaming 
effects, modulation by source rotation, the size of the emit- 
ting region), or even by the properties of space itself (in the 
case of cosmological models). In the past, a frequently used 
timescale has been the duration, which characterizes the 
time interval over which emission is observed above the 
background level. 

The measurement of the duration of high-energy tran- 
sients has historically been difficult. The measured dura- 
ions of identical bursts emitted from identical sources 
located at very different distances could be substantially 
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different, being affected by the signal-to-noise ratio of the 
observation. Comparing duration measurements between 
various experiments can be further complicated by differ- 
ences m the sensitive energy ranges, integration timescales 
and the local background flux (this will be addressed further 
in § 4). 

We will discuss the use of two parameters, T go and T, n 
Kouve hotou et al. 1993), as measures of the characteristic 
timescale over which transient sources emit radiation (i.e. 
the duration ”). In § 2 we provide a definition of T qo ( T . n )' 
an extensive description of how these quantities are mea- 
sured for GRBs observed with BATSE, and a discussion of 
their statistical uncertainties ST go and <5T 50 . In § 3 we will 
explore the systematic effects associated with the measure- 

mC cl ° f Tgo ^ Tso using both simulated simple GRB time 
profiles, and observed burst profiles whose signal-to-noise 
ratios are varied. The methods used to simulate GRBs and 
to vary the signal-to-noise ratio of observed bursts are dis- 
cussed in detail in this section. In § 4 we discuss T qo (T< n ) 
distributions in general, and how they can vary between 
dinerent experiments. We then briefly review the BATSE 
T 90 (T 50 ) results, followed by a discussion of the impact of 
the systematic effects on these BATSE observations. We 
t en review recent T 90 and T 50 measurements made with 
other experiments. Our concluding remarks are presented 

datcc • pendlX A P rovides pertinent details about the 
BATSE instrumentation and data types. Appendix B dis- 
cusses the details of the analytic calculations of T qo (T, n ) for 
the simple time profile shapes used in our simulations 

2. T 90 AND 7^ 0 AS MEASUREMENTS OF DURATION 

2.1. Definition of T 90 and T so 

. 7?°J S de ^ ned as the time it takes to observe 90% of the 
total background-subtracted counts S (01 in a burst, starting 
and ending when 5% and 95% of 5 tol have been observed 8 

sno/ eCt f V fh y ' Similarly ’ T so is defined as the time over which 
50 /o of the counts S lot has been observed, starting and 
ending when 25% and 75% of S to „ respectively, have been 

observed. Clearly, T 90 must always be greater than or equal 
to / 50 . 
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Let us represent a model source count rate history, 
summed over the entire observable energy range, as dS/dt. 
This model is assumed to be free of background. Through- 
out this work we will indicate differential elements with a 
lower-case d, while we will use the Greek symbol A to indi- 
cate discretely binned data, as in the observed count rates 
AC/At. The total source counts observed from a given event 
are then given by 



The function dS/dt is identically zero before the start time 
t s , and after the end time r e , of the transient emission from 
the source. Thus the limits of integration in equation (1) 
may be taken to be from r s to r e . We designate the time at 
which/% of the total source counts have been observed as 
Ty ; it is given by 

f _ % (ds/dt)dt 

100 s tot • (2) 

T go is then given by 

^90 = *95 — T s (3) 

and T so is similarly given by 

^50 = *75 — *25 • (4) 

T go and T 50 are theoretically independent of the burst 
signal-to-noise ratio S/N. The extent to which this is an 
idealization will be discussed extensively in § 3. 

The parameter T go is not equivalent to the duration T d of 
a transient. The 10% of the total emission lying outside of 
the T go interval may be weak and extended (as with time 
profiles that have a fast rise followed by a nearly exponen- 
tial decay), causing a significant difference between the 
values of duration and T go . 

The duration can also differ from T go by the way that the 
two quantities characterize events exhibiting temporal 
structures separated by extended periods with no observ- 
able emission above background. Figure 1 shows a GRB 



Fig. 1. Time history of burst IB 910522. A weak precursor at f % 0 s 
is found to contain less than 5% of the total counts in the event, thus being 
excluded from the T 90 interval. 


observed with BATSE in which there is a single pulse 
(containing less than 5% of the counts in the entire event) at 
the trigger time, followed by a much larger emission 
episode. The T go interval was found to start at ~ r trg + 108 s 
and end at ~r trg + 138 s, giving a T go of ~30 s. However, 
by eye one would estimate the duration 7^ to be ~ 140 s. 

2.2. Calculating T 90 and T so from the BA TSE Data 

The BATSE data used for calculating T go and T 50 are the 
discriminator data from the Large Area Detectors (LADs) 
(the BATSE instrumentation and data types are briefly 
described in Appendix A; for more detailed information see 
Fishman et al. 1989 and Horack 1991). These include the 
DISCLA data, with a time resolution of 1.024 s, the PREB 
and DISCSC data, with a time resolution of 64 ms, and 
sometimes the TTE data, binned to a time resolution of 1 
ms. Each of these data types consists of the count rates 
observed in the four LAD discriminator energy channels 
separately; all data used here are summed (either onboard 
or on the ground) over the triggered detectors. 

For each burst, the data types are concatenated to span 
the time range from ~ 1 15 s prior to the burst trigger time 
T trg to ~r trg + 360 s during the early part of the mission 
(before 1992 December 17), and out to ~z lrg + 550 s after- 
ward. Each time interval in this range is covered by the data 
type offering the finest time resolution available. Thus the 
concatenated data usually consist of DISCLA data from 
~r, rg — 1 15 s to t lrg — 2 s, PREB data from ~r lrg — 2 s to 
*trg> DISCSC data from i (rg to r Irg + 240(550) s, and during 
the early part of the mission the interval from ~ r 1Ig + 240 s 
through r, rg + 360 s is spanned by DISCLA data once 
again. The time range of the concatenated data can be 
extended in either direction using DISCLA data if the tran- 
sient event displays either pretrigger emission, or emission 
that extends significantly later than r lrg + 240(550) s. If the 
event has a short duration (<2 s), or if it contains a rapid 
rise in intensity near / = t (rg , a short time interval near the 
trigger time may be covered by TTE data rather than 
PREB and DISCSC. The start and end times of this interval 
depend on the intensity of the burst (see Appendix A for 
more details), but this interval typically starts no earlier 
than T lrg - 0.18 sand ends no later than ~r trg + 2 s. 

We designate the count rate observed in the ith energy 
channel during theyth time bin as AC(E„ tj)/Atj, where At, 
is the duration of theyth time bin. This indexing convention 
will be used for the remainder of this work. The omission of 
indices will imply that the entire range of that index should 
be considered. Figure 2 a shows an example of a GRB time 
history, made using the concatenated data types, summed 
over all four BATSE energy channels ( ~ 20-2000 keV). 

The background count rates, which we will represent as 
dB{E, t)/dt, in each of the four discriminator channels must 
be modeled before their subtraction from the observed 
count rates. Background models are created by using back- 
ground time intervals during which the source is believed to 
be below the LAD detection threshold, fitting the observed 
data in these intervals to a polynomial of order p , and inter- 
polating the fit across the source intervals. We used values 
of p between two and four; a different index is allowed for 
the background model in each channel. Any number of 
background intervals can be used, though one interval on 
each side of the event is often sufficient. When the burst 
consists of episodes of flux separated by large time intervals 
during which the source appears to be absent, it is some- 
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Fig. 2. — (a) Time profile using concatenated data types, summed over 
the energy range -20-2000 keV. for burst 3B 940217. (b) A plot of the 
integrated counts as a function of time for the burst shown in (a), after the 
subtraction of a polynomial background model from the observed count 
rates. 


times helpful to constrain the model fit by choosing a back- 
ground interval during the quiescent period. It is important 
that the background intervals do not overlap with source 
emission. One should be particularly careful when the time 
history contains weak flux that is changing on a timescale 
similar to that of the orbital variation of the background 
rates. Background modeling is also critical for low-fluence 
events, where fluctuations in background may have fluences 
comparable to the total event fluence. In defining back- 
ground intervals, we have also been careful to avoid time 
intervals during which the count rates contain contribu- 
tions from an active discrete source exhibiting pulsed flux 
(such as Vela X-l) or X-ray flickering (as is often observed 
from the black hole candidate Cyg X-l), or when intensity 
steps due to occultation of a source by the limb of the Earth 
occur. It is also important not to contaminate background 
intervals with transient flux from unrelated GRBs, solar 
flares, or magnetospheric particle events. 

The background model A B(E, t)/At is subtracted from 
the observed count rates AC(£, t)/At to give background- 
subtracted count rates AS(E, t)/At in each channel. The 
integrated source counts S(t) observed during any time 
interval ( t a , t) is given by 


S(f) = I 




AS(tj) 


At , 


(5) 


where AS(r)/Af represents the background-subtracted count 
rates summed over all discriminator channels and_/ is a time 
bin index. We compute the integrated counts as a function 
of time over a user-defined time interval (/i , tj 2 ), where tj 2 
is chosen before the start of the event and t J2 is chosen after 
the end of the event. As an example, the solid curve in 
Figure 2b shows the integrated counts (summed over all 
four energy channels, with tj 2 = — 100 s and tj 2 = 300 s) as 
a function of time for the burst in Figure 2a. Note that the 
data must be complete over the interval (fjj, tj 2 )', S(f) given 


by equation (5) is unknown for all times after any missing 
data. Thus, T 90 (T ?0 ) cannot be calculated for any event 
during which there is a data gap. 

If our background model were perfectly accurate, then 
the solid curve in Figure 2b would remain flat from the first 
point until the first detection of source emission. This time 
interval, immediately preceding the onset of the burst emis- 
sion, will hereafter be referred to as the zero-fluence interval. 
As burst emission is detected, the slope at any point on this 
curve indicates the intensity of the burst. The curve would 
ultimately return to a line with a slope of zero, indicating 
the cessation of emission from the source. The time interval 
following the last detected source emission will hereafter be 
referred to as the total-fluence interval. 

Theoretically, negative slopes in Figure 2b are not pos- 
sible; however, in practice negative slopes result from dis- 
crepancies between our background model and the true 
background count rates. These discrepancies may occur 
over long timescales due to the observed background rates 
varying faster than variations allowed by our low-order 
polynomial models (for example, when the spacecraft is in 
the proximity of the South Atlantic Anomaly), or they may 
occur over short timescales, due to Poisson fluctuations in 
the counting rates. 

The start and end times of the zero-fluence and total- 
fluence intervals are manually defined; these intervals 
should not overlap with any source emission. The zero- 
fluence level L 2 , and the total-fluence level L ( , represent the 
value of the integrated counts S(t) immediately prior to the 
onset of the source emission, and immediately after the cess- 
ation of the source emission, respectively. L z and L, are 
found by calculating the mean value of S(t) during the user- 
defined zero-fluence and total-fluence intervals, respectively. 
The first time at which the integrated counts S(f) exceeds L z 
(L,) is considered the start (end) time of burst emission. The 
spreads of the data about the mean in the zero-fluence and 
total-fluence intervals are characterized by their variances 
{dL z ) 2 and (dL,) 2 , respectively. 

The total number of counts AL observed due to the GRB 
is then 

AL = L, - L z , (6) 


and the/% fluence level, S f , is given by 


/ 


s '- L - + ™ aL 


( 7 ) 


Sf represents the value of the integrated counts S(t) when 
/% of the total source counts have been detected. The time 
at which the integrated counts reaches the value S s can then 
be found and designated as z f . The true value of z f can fall 
anywhere within a data bin. When the appropriate data bin 
in which z f is contained has been found, then the value 
assigned to Zf is the start time of that bin if/is either 5 or 25, 
or the end time of that bin if/is either 75 or 95. The values 
of T 90 (T 50 ) are then found from equations (3) and (4). 


2.3. Statistical Uncertainties in T 90 and T 50 

The uncertainties in S f , (dS f \ ot , consist of contributions 
from two components: ( dS f ) czlV due to the uncertainty in the 
integrated counts S(f) at any time t caused by detector 
counting statistics, and (dS f ) fWc , due to the uncertainty in 
the levels L z and L, caused by statistical fluctuations about 
the smooth background model. The uncertainty (dS f ), ot is 
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given by 

(ds f \ ot = . (8) 

Assuming no uncertainties in our background model 
AB(E, t)/At, the uncertainties (SS/dt)j in the background- 
subtracted count rates A S(E, t)/At and the uncertainties 
(5C/St)j in the observed count rates AC(E, t)/At are equiva- 
lent. If we represent the observed count rates in time bin j, 
summed over the four energy channels, as AC(tj)/Atj, the 
uncertainty (dS f ) cnt is then given by 



j = to 


AC(t f ) 

- At i . 


A tj 


( 9 ) 


where r 0 is the time at which the integrated counts exceed 
the level L z . 

From equations (6) and (7), the 5% fluence level S 5 can be 
written as 


S 5 


19L. + L, 
20 


( 10 ) 


The uncertainty (</S 5 ) flllc in S 5 is given by 


(<fS 5 ) f |uc — 



( 11 ) 


where dL z and dL, are the variances associated with L. and 
L„ and the partial derivatives are found analytically from 
equation (10) to be dS s /dL z = 0.95 and dSJSL, = 0.05. 
Equation (11) becomes 

(dS 5 ) nuc = v/(0.95) 2 dD z + (0.05) 2 dL} . (12) 

Similarly, we find 


Sqe — 


L z + 19 L, 
20 


and 


(^ 95 ) fluc = y/ (0.05) 2 dL 2 + (0.95) 2 rfL 2 


(13) 

(14) 


The uncertainties (dS 5 ) lol and (dS 95 ) lol can now be found 
using equations (8) and (9). 

The times and t /+ are the times at which S f - 
{dSf ) tot an <i Sf + ( dS f ) tox counts, respectively, have been 
detected. We define Ax f as 

Ax s = T f+ - Ty_ . (15) 

The statistical uncertainty <5T 90 in T 90 is now calculated by 

&T 90 = \Z(At 5 ) 2 + (At 95 ) 2 . (16) 

It is through the quantities At 5 and At 95 that the depen- 
dence of ST 90 on the detailed shape of the burst time profile 
is introduced. The error 5T S0 in T so is computed similarly. 

Note that the minimum value of Ar f that can be achieved 
is t ICS , where r res is the time resolution of the discrete data at 
the time x f . Assuming that the values of x f fall within data 
bins of equal time resolution t res , the minimum value of the 
uncertainties ST 90 and <5T S0 given by equation (16) is 

tre S (2) 1/2 . 

The method used to calculate the statistical uncertainties 
ST 90 and 5T S0 in the BATSE IB (Fishman et al. 1994), 2B 
(Meegan et al. 1994), and 3B catalogs (Meegan et al. 1995) 
did not include the contribution of ( dS f ) cnt to (dS f ) m . This 
effect is insignificant for short events, where the sum in 


equation (9) contains a small number of terms. After recal- 
culating the uncertainties for a number of bursts, we esti- 
mate that the statistical uncertainties ST go and ST 50 for 
long-duration bursts in these three BATSE catalogs should 
be larger by a factor of ~ 1 .5. 

2.4. Comparison ofT 90 and T 50 

Some believe that T S9 is a better measure of burst dura- 
tion than T 90 , because it is less susceptible to weak flux in 
the wings of the burst (Kargatis et al. 1994). This statement 
is correct only under the assumption that weak flux prefer- 
entially occurs near the start and end of the burst, and 
avoids the intermediate region of the burst where t 25 and 
t 75 are likely to be found. Both T qo and T 50 are susceptible 
to weak flux, and which is affected more depends entirely on 
the detailed profile shape. Similarly, T 90 is not necessarily 
more sensitive than T so to precursor emission separated 
from the main burst emission by an extended background 
interval, unless one assumes a priori that the flux is more 
intense at t = r 25 than at t = t 5 . By examination of Figure 
2b, it is apparent that if the first 5% of the counts observed 
from the burst were removed, the resulting effect on the 
values of Tj is dependent on the shape of the time profile. 

To investigate this issue further, we examined the dis- 
tribution of the relative uncertainties <5 T 90 and ST 50 in the 
BATSE 2B catalog, shown in Figures 3a and 3b, respec- 
tively. The two distributions are very similar, thus not 
lending support to a preference of one measure over the 
other. The average percentage uncertainties are <<5T 90 > = 
19 + 2.2 (36)% and < [ST so y = 16 + 1.5 (24)%, where the 
first error indicates our uncertainty in the mean and the 
error in parentheses is the sample standard deviation and 
characterizes the width of the distribution. The majority of 
the large percentage uncertainties occur at small values of 
T 90 (Tso), and are the result of the finite time resolution of 
the binned data. Not all BATSE measurements of small T 90 
(T 50 ) have large percentage uncertainties. Recall that 
BATSE observations of bursts often include TTE data (see 
Appendix A) which are used at a time resolution of 1 ms. 
For those bursts covered entirely by TTE rather than the 
coarser DISCSC, the minimum uncertainty in T 90 (T 50 ) 
(given by eq. [16]) drops from ~0.091sto ~ 1.4 ms. 


2.5. Other Quantities 

Hardness ratios, which roughly characterize the source 
spectrum, can be obtained from the total counts detected in 
each energy channel during the T 90 (T so ) interval (with the 
interval determined as just described from the rates 
summed over all energy channels). These hardness ratios 
were calculated for the BATSE data and were used to 
confirm that on average short bursts have harder spectra 
than longer bursts (Kouveliotou et al. 1993). 

Two other quantities which may be found are T r and T dcy , 
which characterize the time required for the event to rise to 
peak intensity and to decay from peak intensity, respec- 
tively. The time r pk is the time at which the background- 
subracted count rate, summed over all energy channels, is 
maximal. The time T r to rise from the start of the T go inter- 
val to peak intensity is then given by 

T r = fpk-is- (17) 

Similarly, the time 7^ cy to decay from the peak intensity to 
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Fig. 3 a 

Fig. 3.— (a) The distribution of the percentage uncertainties in T w for 
uncertainties in T 50 for the same sample. 


the end of the T g0 interval is given by 

T acy = t„ - T pk . (18) 

These quantities are different from the usual rise and decay 
times which measure the times to rise from, and decay back 
to, the background count rate. Notice that it is possible for 
either T r or T dcy (though not both) to be negative, indicating 
that the peak intensity of the burst occurred cither before t 5 
or after t 95 , respectively. 

One can define a parameter A pk , which characterizes the 
position of the peak intensity within an event, as 


If one accounts for negative values of T r by setting them 
equal to zero, the value of A pk is constrained to fall within 
the range 0 < A pk < 1. This parameter provides a sense of 
the asymmetry of the overall time history. If the peak inten- 
sity of the event is equally likely to occur at any time within 
the event, we would expect an average value (A pk ) ~ 0*5, 
with an associated uncertainty that is due to a finite sample 
size. 


3. SYSTEMATIC EFFECTS ON T 90 AND T 50 

In this section we explore some likely sources of system- 
atic errors using simulated time profiles, and observed time 
profiles whose signal-to-noise ratios are varied. The simu- 
lated time profiles can be represented by simple formulae, 
and their values of T 90 (T 50 ) can be calculated analytically. 
We verify that the code used to calculate T 90 (T 50 ) works 
properly, by applying it to simulated noise-free burst pro- 
files. We examine the effects of source intensity and back- 
ground modeling on T 90 (T 50 ) for simulated bursts with 
statistical noise. We then systematically reduce the signal- 
to-noise ratio of several observed bursts, and examine the 
dependence of T 90 (T 50 ) on the detailed shape of the time 
profile as well as its intensity. 

3.1. Simulation of Simple Time Profiles 

The simulated data consisted of count rates AC(£, t)J 
At as a function of time, in each of the four BATSE discrimi- 



Fig. 3 b 

the bursts in the BATSE 2B catalog, (b) The distribution of the percentage 


nator energy channels, with a time resolution identical to 
that of the concatenated data described in § 2.2. 

We first created background model count rates A B(E, t) 
/At in each of the four energy channels, by randomly selec- 
ting a set of polynomial coefficients from those used to 
model background in the T go calculations of the 2B catalog. 
These coefficients were then used for the analyses discussed 
in §§ 3.1.1 and 3.1.2 (the effect of using different background 
models is considered in § 3.1.3). The only variations in this 
background model are long-term and are due to the orbital 
motion of the CGRO about the Earth. Statistical fluctua- 
tions in the count rates have not yet been added. 

To simulate the source count rates A S(E, t)/At we first 
had to determine how to divide the total counts observed at 
any given time between the four discriminator energy chan- 
nels. To accomplish this, we used a 2B catalog database 
consisting of the total counts observed in each energy 
channel during the T 50 time interval. We neglect in this 
work the systematic effects on T g0 (T 50 ) that may be intro- 
duced by varying spectral shapes; therefore, one count spec- 
trum (shown in Fig. 4) was randomly selected from the 2B 
catalog and it was used for all of the burst time profile 
simulations. The shape of the count spectrum remains con- 
stant throughout each simulated burst; no spectral evolu- 
tion is allowed. 

We simulated three simple burst profile shapes: triangu- 
lar, rectangular, and one exhibiting a fast rise followed by a 
nearly exponential decay (hereafter referred to as a FRED). 
These three forms were chosen because each could be 
described in a straightforward manner by a simple analytic 
expression, which allowed for their values of T go and T 50 to 
be calculated analytically. 

The analytic model (dS/dt ) f of a triangular time history in 
each energy channel may be written as 



0 , 

m, t + b, , 

n, t + d; , 
0 , 


r < t s , 

T s < t < Tpk , 
Tpk < t < T,. , 

t > z e , 


( 20 ) 
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Fig. 4. — Count spectrum of BATSE burst IB 911129, integrated over 
the T 50 interval; this spectrum was used in the simulations of burst time 
profiles. All values are normalized to the value observed in the first energy 
channel (20-50 keV). 
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where z s is the start time of the source emission, i pk is the 
time of the peak emission, and x e is the end time of the 
source emission. Each of these three time constants are 
input in seconds relative to t = 0, and are independent of 
energy channel, for each simulation. The slopes m k and of 
the burst rise and decay are given by 


nii = 



and 


n i = 



( 21 ) 


where P t is the peak count rate in energy channel i. The 
constants b { and d, are given by 


hi = -x s mi and d { = -x e n i . (22) 

Note that the energy dependence of equation (20) is entirely 
through the peak count rate P, in each energy channel. 

The analytic model of a simulated rectangular time 
history in energy chanel i may be written as 

(dS\ _ f 1 3 • 

U/, ] Pi> T * - f < T - ( 23 > 

[0 , t > x e , 


where x s and x e have the same meanings as before. As with 
the triangular profile, the energy dependence is provided 
through the peak count rate P, in energy channel 

Our FRED profile is an idealization of a subclass of 
bursts that are actually observed in the data set 
(Kouveliotou et al. 1991), The analytic expression for the 
model used to simulate a FRED in energy channel i is given 
by 


dS\ = f° . t < t, , 

dtji {P i e~’“ , t> t s , 


(24) 


where t s is the start time of the emission, x d is the \/e decay 
time, and k = x d *. As before, P ( provides the energy depen- 
dence. Our simulated FREDs differ from those observed in 
the BATSE data set in that the observed bursts do not 
usually rise so rapidly to their peak intensity, and the decay 
of the observed bursts is seldom precisely exponential. 


The final parameter of equations (20) through (24) that 
needs to be determined is the peak count rate P, in energy 
channel i. The quantity used to characterize the burst inten- 
sity is the signal-to-noise ratio S/N at the peak of the time 
history summed over the energy channels 2 and 3. Using the 
background model we calculate the background counts 
(# P k)2 + 3> summed over channels 2 and 3, observed at the 
peak, and the “noise” (t 2 + 3 at the peak time, given by 
[(£ p k)2 + 3] I/2 - The source counts (N pk ) 2 + 3 at the peak, 
summed over channels 2 and 3, corresponding to an input 
S/N is then given by 

(Wpk)2 + 3 =(S/N)<7 2+3 . (25) 

These (N pk ) 2 + 3 counts must be divided between channels 2 
and 3 according to the previously selected count spectrum. 
Using the counts in these two channels and the assumed 
count spectrum, the peak count rate P, in energy channel i is 
straightforwardly determined. 

With the peak count rate P, determined in each energy 
channel, and each of the appropriate time constants defined, 
the source count rates A S(£, t)/At in each channel are then 
calculated according to equation (20), equation (23), or 
equation (24). The total observed count rates AC'(£, t)/At 
are then determined by summing together the background 
model count rates and the source count rates, 

AC'(£, t) AP(£, t) A S(£, t) 

At ~ ~~aT~ + ’ (26) 

where the primed notation indicates that no statistical noise 
has been added to these rates. Therefore at this point the 
input signal-to-noise ratio S/N does not yet have its usual 
meaning, and has only been used as a measure of the burst 
intensity above a smooth background model. However, 
these rates can be very useful in testing the T 90 (P 50 ) soft- 
ware (see § 3.1.1). 

The final step in the simulation of the simple burst pro- 
files AC(£, f)/Af is to add statistical noise to the count rates 
AC(£, t)/At in each energy channel. The number of counts 
AC in the 7th time bin with noise added is given by 

ACj = A Cj + Aipj , (27) 

where Ai/^ represents the statistical noise and is randomly 
generated from a Gaussian distribution with a mean // = 
A Cj. The S/N now has its traditional meaning, and is more 
than just a measure of the burst intensity. The total count 
rates AC(£, t)/At in each energy channel are then given by 

AC(£, 0 A £(£, t) ( AS(£, t) t A*(£, t) 

At ~ At (28) 

After each simulation we checked whether the burst 
would have set the BATSE burst trigger on at least one of 
the three trigger timescales, under the assumption that the 
burst was observed in two LADs, with equal intensity in 
each detector (see Appendix A). 

Figure 5 shows three sample simulated bursts with a 
background model and statistical noise added. Figure 5a is 
an example of a triangular burst, with S/N = 50 o at the 
peak, and duration T d = x e — x s = 30 s. Figure 5b shows an 
example of a rectangular burst, with S/N = 5.5 o and a 
duration T d = 1.5 s. Figure 5c shows an example of a 
FRED, with S/N = 1 5 a and a decay time x d = 70 s. 
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Time (seconds since trigger) 



Fie.. 5 —Examples of simulated burst time histories, with a background 
model and statistical noise added, (a) Triangular burst, {h) Rectangular 
burst, (t ) Burst with a fast-rise and an exponential decay (FRED). 


3.1.1. Using Simulated Time Profiles to Test the T go Code 

We used the procedure described in §§ 2.2 and 2.3 to 
calculate T 90 (T so ) for a series of simulated bursts which 
were free of statistical noise. We then analytically calculated 
T 9Q (T so ) for each of these simulated bursts (these calcu- 
lations are described in Appendix B). Nine simulations were 
run for each of three simple time profile shapes: triangular, 
rectangular, and FRED. The nine simulations included 
bursts with durations T d of 1.5, 30, and 150 s (or decay times 
z d of 1.5, 30, and 70 s in the case of FRED profiles). For each 
of these profiles we assumed three different intensities (S/N): 
weak (5.5 cr), moderate (15 er), and strong (50 o). The signal- 
to-noise ratio S/N is still defined by equation (25); though 
noise has not been added to these profiles, o 2 + 3 can still be 
calculated. These signal-to-noise ratios do not have their 
traditional meaning; they are simply measures of the burst 
peak intensity. 

The results of these simulations are given in Table 1, in 
the form of the largest relative error (E R ) as compared to the 
analytic values. Statistical uncertainties on the measured 
values cannot be calculated for these simulations due to the 
lack of statistical noise in the time profiles. 

We see from Table 1 that significant errors occur for 
profiles with very short durations. These errors occur 
because the BATSE time histories consist of discretely 
sampled data while the analytic calculations are done using 
continuous functions. The measured values differ from the 
analytic values by less than one data bin. This is why the 
percentage errors increase as T d decreases. This also 
explains why the percentage uncertainties in T 50 are larger 
than those in T 90 ; T so is always smaller than 7^ 0 and there- 
fore is more susceptible to binning effects. 

For each triangular burst of the same duration, the three 
measured values of T 90 were identical, independent of the 
burst intensity. The same was true for the measured values 
of T 50 . This was also true for the rectangular bursts. In each 
case the measured values of T 90 and T 50 for the triangular 
bursts were larger than the analytic values. The three mea- 
sured values of T 90 for the rectangular bursts with T d = 1.5 s 
were lower than the analytic values; all other values of T 90 
and T 50 for the rectangular bursts were larger than the 
analytic values. For the short FREDs, all values of T 90 and 
T 50 were larger than the analytic values, but by less than 5% 
and 2%, respectively. For these shorter FREDs, the values 
measured for T 90 and T 50 were independent of the burst 
intensity. 

The results from the simulations of FREDs with longer 
decay times are indicative of the difficulty of making the 
measurements when the burst flux varies on timescales 
similar to those over which the background model is 
varying. Figures 6 a and 6b show the measured values of T 90 
and Tso, respectively, as a function of the simulated burst 
intensity for FREDs with z d = 30 s. The dashed lines indi- 
cate the known analytic values. It is clear that T 90 and T 50 
show no systematic dependence on the burst intensity; 
however, the measured values of T 90 and T 50 are systemati- 
cally lower than the analytic values by ~ 2% and ~1%, 
respectively. 

The systematic errors are especially pronounced in the 
measured T 90 and T so values of the simulated FRED pro- 
files with z d of 70 s (see Figs. 6c and 6 d). We attribute these 
systematic errors to inadequate modeling of the back- 
ground rates, largely due to the fact that the tail of the decay 
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TABLE 1 

Simulation Results for Noise- Free- Profiles 



Triangular 


Rectangular 


FRED 



^90 

^50 

^90 

Tso 


^90 

To 

T d 

Er 

E „ 


Er 



Er 

( s ) 

(%) 

(%) 

(%) 

(%) 

( s ) 

(%) 

(%) 

1.5 ... 

<7 

<17 

<1 

<12 

1.5 

<5 

<2 

30 ... . 

<1 

<2 

<1 

<1 

30 

<2 

<1 

150 ... 

<1 

<2 

<1 

<1 

70 

<18 

<10 


in intensity of the FRED is varying on the same timescale as 
are the background rates, and it is then very difficult to 
separate the two (even in this noise-free case). However, we 
do not find any systematic dependence of either T go or T so 
on the burst intensity. The systematic effects of background 
modeling and the detailed background model on T 90 (T so ) 
will be addressed in § 3. 1 .3. 

We would like to point out that the simulated triangular 
burst with 7^ = 150 s and S/N = 5.5 would not have trig- 
gered the BATSE instrument. Therefore it was not included 
in Table 1. It will be seen in § 3.1.2 that when we simulated 
the same burst with statistical noise added, we found that 
the instrument would have triggered. 

Using noise-free GRB simulations, we find that the 
primary sources of error in these measurements are the 
finite width of the discretely sampled data and the uncer- 
tainty in background modeling; the first error is unavoid- 
able, and the second is examined in greater detail in § 3.1 .3. 



3 . 1 . 2 . Dependence ofT go on Signal-to-Noise Ratio 

We use the standard BATSE software to calculate T 90 
(T»o) for a series of simulated bursts which now contain 
statistical noise. Fifteen simulations were run for each of the 
three simple time profile shapes. The 15 simulations 
included bursts with durations T d of 1.5, 30, and 150 s (or 
decay times x d of 1.5, 30, and 70 s for FRED profiles). For 
each duration, profiles of five values of S/N were simulated: 
5.5, 10, 15, 30, and 50 a. In these simulations, the signal-to- 
noise ratio has its traditional meaning. The results of these 
simulations are given in Table 2, in the form of the largest 
relative error (£„) and the largest absolute deviation E A 
from the analytic value. We note that for each burst profile, 
the size of the uncertainties ST go ( 5T 50 ) was inversely pro- 
portional to the burst S/N. 

For the triangular and rectangular bursts, measured 
values of T go (T so ) that were less than the analytic values 




1 10 
S/N (a) 


10 

S/N (a) 


an^y^icall/calcluat^d^^u^^f T 9 /and S T 5 ^ 0 ^fl) r" ^=^0 sjc) ^ UrSt S/N f ° r S,mU ' a,ed n ° ,Se - free FRED P rofiles ' ind.cate the 
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were observed, but on average the measured values were 
larger. In no case did we find evidence for any systematic 
dependence of either T 90 or T 50 on the burst S/N. We find 
no evidence in the FRED case of t, = 1.5 s for any depen- 
dence of either T go or T s0 on the burst S/N. 

We find that the largest errors are associated with the 
longer duration FRED profiles (see Figs. 7 a and lb). For 
these FREDs with z d = 30 s, the errors in T go and T 50 were 
les than 20% and 10%, respectively. The measured values 
are on average underestimated ; however they are consistent 
with no systematic dependence of T go or T 50 on the burst 

S/N. . 

The systematic errors are pronounced in the measured 
T g0 and T 50 values of the simulated FRED profiles with 
T = 70 s (see Figs. 7c and Id). A systematic dependence of 
T g o and T 50 on the burst S/N is apparent, the bursts with 
smaller S/N having larger deviations of T go and T 50 from 
their analytic values. We will show in § 3.1.3 that this sys- 
tematic dependence on burst S/N can be attributed to the 
inadequate modeling of the background rates. 

If T go (T 50 ) were systematically dependent on the burst 
S/N, then one would expect a tendency for the measured 
values to be consistently overestimated or underestimated 
in the weak case (S/N = 5.5 a) for all profile shapes. The 
values of T go and T 50 for the weak triangular bursts were 
overestimated. There was no tendency for either underesti- 
mation or overestimation of 790 (Tso) rec- 

tangular bursts. T g 0 (T so ) were underestimated in the cases 
of the weak FRED bursts, except for the case of z d = 1.5 s; 
however, we attribute this underestimation to incorrect 
background modeling. Thus, we conclude that any system- 
atic variations of T g0 (T 50 ) with burst S/N depend on the 
shape of the burst time profile. 

We point out that in the case of one simulated triangular 
burst, that with T d of 1.5 s and a S/N of 30, the uncertainty 
6T g0 in the measured value of T go was large, given by T g0 = 
1.1 + 9.7. This is due to the chance superposition of positive 
statistical fluctuations in the count rates in energy channels 
1, 2, and 3. When the counts in these fluctuations were 
summed, they resulted in a large value of t 95 + (defined in 
§ 2.3), thus a large value of At 95 given by equation (15), and 
consequently a large uncertainty ST go . This is an indication 
that the value of the total fluence level L, (see § 2.2) was not 
well determined. This illustrates the fact that the total 
counts in background fluctuations can be significant, com- 
pared to the total counts observed in short bursts, and can 


result in large statistical errors on measurements of T go 
(T 50 ) for these bursts. 

We also notice that unlike the simulated noise-free burst 
with S/N of 5.5 and T d of 150 s, the same simulated burst 
profile with statistical noise added would have triggered the 
BATSE instrument. 

3.1.3. Dependence of T go on Background Effects 

The effects of inadequate background modeling on T g0 
(T 50 ) measurements are explored in this section. We will 
first investigate the effect of defining postburst background 
intervals during which there is still burst emission. This is of 
particular concern in events where there is low-level flux, 
varying on timescales similar to those of orbital variations 
in the background rates. We will then address the depen- 
dence of T go (T 50 ) on the shape of the background. 

Two simulated burst profiles were used to examine the 
effects of defining the postevent background interval before 
the end of the event. The burst profiles were FREDs, with a 
decay time z d of 70 s and signal-to-noise ratios S/N of 50 a 
and 5.5 a. The same background model, a fourth degree 
polynomial, was used for each of the two bursts. The count 
spectrum used to simulate the burst count rates in each 
channel is shown in Figure 4. 

To characterize how near the postevent background 
interval is to the burst, we define the parameter 

x = ffcl ~ t .21 ( ( 29 ) 

T 90 

where t bl is the start time (relative to t = 0) of the postevent 
background interval, and the analytic value of T go given by 
equation (40) in Appendix B is used. The parameter x mea- 
sures how much later, in units of T 90 , the start time of the 
postburst background interval is than the end time of the 
T go interval. Low values of x may lead to problems of 
including source flux in the background fits, and high values 
of x may not adequately constrain the background model 
during the time of the burst (particularly, as we will see, 
when high-order polynomials are used to characterize the 
background). 

We first considered the 50 a FRED with x d = 70 s. By 
varying t b j, we measured T go and T so for nine different 
values of x, approximately equally spaced between ~0.0 
and 0.8. Though background with a polynomial index of 4 
was used during the simulation of the burst, the back- 


TABLE 2 


Simulation Rlsults for Time Profiles with Noise 



Triangular 


Rectangular 


FRED 


T d 


^50 

eje a 

^90 
Er/E A 

To 

Er/E a 


^90 

Er!E a 

To 

E r IE a 

(s) 


(%/<r) 

(%/ff) 

(%/a) 

(s) 

i%/o) 

(%/< d 

1.5 

< 12(45)71 
<6/1 

<12/1 

<5/1 

<11/1 

1.5 

<35/l(2) b 

<20/1(2/ 

^0 

<2/1 

<2/1 

<2/l(2) d 

30 

< 20/* 

< 10/ c 

150 

<2/2(4/ 

<3/2 

<1/2 

< 1/2 

70 

< 20/ c 

<10/ e 


• The S/N = 5.5 a case had £ K * 45%, but E A < l a. 
b The S/N = 10 a case had E A < 2 a and £„ a 35%. 
c The S/N = 15 a case had E A <2a and £„ a 15%. 

“ The S/N = 30 a case had E„<2a and £„ a 1.5%. 
c The measured values are more than 5 a from analytic value. 
f The S/N = 10 a case had E A < 4 <r and £„ a 2%. 



No. 2, 1996 


DURATION OF GAMMA-RAY BURSTS 


579 




220 

r 

i ii r ni| 1 — i i mu 

200 

- 

O 

180 



160 

- 

- 

140 

- 

(c) - 

120 

l_ 

1 1 1 1 ml 1 1 » 1 1 Hi 



7 - ^ ^ ^he measure d values of T 90 and / 50 as a function of burst S/N for simulated FRED time profiles to which statistical noise has been added. 
The dashed lines indicate the analytically calculated values of T 90 and T 50 , {a) and (c) as in Fig. 6. 


ground was modeled with a polynomial index of 2 during 
each calculation of T 90 (it is shown later in this section that 
low-order polynomials should be used to model back- 
ground data that are separated by large intervals of source 
emission). Figures 8a and 8h show the measured values of 
T 90 and T 50 as a function of x. The measured values are 
systematically smaller than the analytic values (which are 
indicated by the dashed lines). The errors in T 90 and T 50 
were never larger than 30% and 20%, respectively, and in 
most cases were smaller than these values by a factor of ~2. 
There is also a clear dependence of the measured values of 
T go and T 50 on x. The nearer the start of the background 
intervals to the source emission, the greater the chance of 
including weak flux in the background model, and the 
smaller the measured value of T go and T 50 . 

We similarly measured T 90 (T 50 ) of the 5.5 a FRED with 
x d = 70 s, using the same nine values of x(t M ) as before. A 
polynomial index of 2 was again used for the background 
model. The majority of the measured values (Figs. 8c and 
8 d) were smaller than the analytic values. The errors in T go 
and T 50 were never larger than 33% and 20%, respectively. 
In spite of larger fluctuations, there is again a clear depen- 
dence of T go (T so ) on x; T go and T 50 become smaller the 
nearer the background interval gets to the event. 

To investigate the effects of using a more flexible back- 
ground model to represent the same background data, we 
reanalyzed both FREDs, using a polynomial of index 4. The 
errors are in general worse than for those when the poly- 
nomial index of 2 was used. For the 50 a FRED (see Figs. 
9 a and 9b ), T go and T so were less than the analytic values by 


up to 35% and 25%, respectively. The systematic trend for 
lower values of T go and T 50 to be derived when lower values 
of x are used is still apparent. For the 5.5 a FRED (see Figs. 
9c and 9 d), T go and T 50 were less than the analytic values by 
up to 55% and 45%, respectively. The systematic depen- 
dence of T go and T 50 on x that is clear in Figures 8, 9a, and 
9b is not seen in Figures 9c and 9 d. This is because at larger 
values of x, t bl is so near the end of the available data that 
there is sometimes not enough data in the postburst back- 
ground interval to adequately constrain high-order poly- 
nomial fits. The resulting fits allow too much variation in 
the interpolated model during the source interval, thus 
causing improper background subtraction; this effect is 
clearly more pronounced for weak, long-duration events. 
We stress that this is caused by the limited amount of simu- 
lated data, which ends ~360 s after the trigger time. This 
does in fact mimic the observed data from early in the 
mission, but during the later part of the mission these data 
extended out to ^550 s, minimizing this effect for all bursts 
other than those of extremely large duration. Nonetheless, a 
lower order polynomial (polynomial index of 2 or 3) should 
be used to fit background data that are separated by large 
periods of source emission. 

We next used simulated bursts to investigate the effects of 
various background shapes on T 9 o (T 5 oI We randomly 
selected 10 background models, whose polynomial indices 
varied between 2 and 4, from those that were used in the 
BATSE 2B catalog (see discussion of these models at the 
beginning of § 3.1). We added the same simulated burst to 
each model. Any discrepancies between the measured and 
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Fig. 8. — (a)-(</) The measured values of T 90 and T 50 as a function of x (defined in the text by eq. [29]) for simulated FRED profiles with r d = 70 s. A 
polynomial of index 4 was used to simulate the background data, but an index of 2 was used to model background during each calculation. The dashed lines 
indicate the analytically calculated values of T 90 and T 50 . (a) S/N = 50 cr; (c) S/N = 5.5 a. 
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p IG 9 — (a)-(d) The measured values of T 90 and T 50 as a function of x (defined in the text by eq. [29]) for simulated FRED profiles with r d = 70 s. In 
contrast to the data shown in Fig. 8, a polynomial of index 4 was used to simulate the background data, as well to model background during each calculation. 
The dashed lines indicate the analytically calculated values of T 90 and T 50 . (a) and (c) as in Fig. 8. 
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analytic values would now be a result of the effects of differ- 
ent background shapes (combined, of course, with statistical 
fluctuations inherent to the measurement process itself). We 
repeated this procedure with four simulated profiles: a 
strong and weak triangular burst (both with T d = 100 s), 
and a strong and weak FRED (both with i d = 30 s). The 
same 10 background models were used for each of the four 
simulated time profiles. When measuring T 90 for these 
simulations, we used polynomials of index 2 to fit the back- 
ground rates, independent of the index used in the simula- 
tion of the rates. This was done to ensure that the 
background model was sufficiently constrained by the data. 
As a result, the scatter of the measured values about the 
analytic values was slightly greater for those profiles which 
had background rates simulated with a 4th degree poly- 
nomial. 

Figures 10a and 10b show the measured values of T 90 and 
T 50 , respectively, for the strong triangular burst with S/N = 
50 a. The x-axis values in Figure 10 have no physical 
meaning, and serve only to identify distinct measurements. 
The measured values of T go (T so ) are scattered about the 
analytic values (indicated by the dashed lines in Fig. 10), 
with no tendency for either underestimation or over- 
estimation. The errors in T go and T 50 were less than 1%. 
The errors (see Figs. 10c and 10 d) associated with the 
weaker triangular burst, with S/N = 5.5 cr, were less than 
8%. The scatter of the measured values about the analytic 
values is larger for the weaker burst, but we do not find any 
tendency for either underestimation or overestimation. 

The errors in T 90 and T so (see Figs. 1 la and 1 lfi) for the 
strong FRED, with S/N = 50 ct, were less than 4%, except 


for one measurement of T go which had an error of ~ 10%. 
There is a slight trend for the measured values to be under- 
estimated; however, this can be attributed to the previously 
discussed effects of large values of x, and not to the varying 
background shapes. The errors in T 90 and T so (Figs. 11c 
and 1 Id) for the weaker FRED, with S/N = 5.5 a, were less 
than 30% and 20%, respectively. Again, the scatter of the 
measured values about the analytic values is larger for the 
weaker burst. We thus find no significant dependence of the 
measured values of T 90 (T 50 ) on the shape of the back- 
ground. 

3.2. Reduction of the Signal-to- Noise Ratio of Observed 
Gamma-Ray Burst Profiles 

To gain insight into the systematic effects associated with 
the detailed burst profile shapes, we chose several bursts 
observed with BATSE, systematically reduce their intensity, 
and in each case measure their T 90 and T 50 . We will first 
describe the details of the signal-to-noise reduction of the 
observed bursts; our results will be discussed in § 3.2.1. 

The time histories of each observed burst consist of LAD 
count rates AC(E, t)/At in the four BATSE discriminator 
energy channels, with a time resolution as described in § 2.2. 
Background data before and after the burst are fitted to a 
polynomial function of order between 2 and 4, and the 
background model AB(E, t)/At is created by interpolating 
the fit across the burst interval. Each energy channel is 
modeled independently. This background model is saved 
and is used for all reductions of the S/N of this burst. 

The source count rates AS(£, t)/At are found by subtrac- 
ting the background model rates from the observed count 



in k":J° V'r' measared values of T 90 and 7~ 50 for two different simulated triangular profiles, both with a duration T. = 100 s using the same set of 
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analytically calculated values of T 90 and T so . (a) and (0 as in Fig. 8. 


rates in each channel, 

A S(E, t) = AC(E, t) _ A BjE, t) (3Q) 

Af _ At At 

The counts AC(E ; , t y ), AS(£„ t } ), and A B(E„ tj) in each time 
bin are calculated. 

We use the count rates AS(E, t)/At as a template for the 
source rates, thus making the assumption that these rates 
are free of any statistical fluctuations. Note that we use the 
same template to create the time histories over the entire 
range of intensities, otherwise we would magnify the sta- 
tistical fluctuations that are inherent to AC(E, f)/Af. 

The time t pk of the peak count rate was determined using 
the data summed over energy channels 2 and 3. Assuming 
that the background counts are Poissonian, the statistical 
noise <r pk at the peak time is given by 

<r pk = v / AB^> 0D 

where AB pk is the number of background counts observed 
during the data interval containing the peak rate. The 
signal-to-noise ratio of the peak, (S/N) pk , is then given by 

(S/N) pk = ^ , (32) 

^pk 

where AS pk is the number of counts observed above back- 
ground during the data interval containing the peak rate. 

The signal-to-noise to which the burst is to be reduced is 
designated as (S/N) r . The reduced source counts A S r (E, t ) 
are then found from 

AS r (E„ tj) = £AS(Ej, tj) , (33) 


where C is defined as 


(S/N) , 

(S/N). 


Note that this process involves reducing the source counts 
in each bin by the same multiplicative factor, and does not 
involve subtracting a constant from each bin. Thus the ratio 
between the source counts in any two time bins will remain 


constant. 

The reduced total count rates AC r (E, t)/At are produced 
using the same method described in § 3.1 (see eqs. [26] 
through [28]). The reduced source rates are first combined 
with the background model, then an appropriate level of 
statistical noise is added. After each simulation we con- 
firmed that BATSE would have triggered at some time on 
at least one of the three timescales. 


3 . 2 . 1 . Dependence of T 90 on Signal-to- Noise Ratiofor V'arious 
Time Profile Shapes 

We applied the procedure described in § 3.2 to a group of 
12 GRBs, (shown in Figs. 12 and 13, summed over channels 
1 through 4) observed with BATSE. We tried to include a 
variety of profile shapes in this group. The values of T go 
(T 50 ) cannot be calculated analytically for these observed 
bursts. Instead of quoting the relative errors and absolute 
errors as in § 3.1.1 and 3.1.2, we will examine relative and 
absolute differences ( D R and Df), with respect to the values 
measured for the originally observed unreduced burst. 
Table 3 lists the values of peak S/N to which each burst was 
reduced, as well as the maximum values of D R and D A for 
both T 90 and T 50 . 
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Fig. 12.— Time profiles of six gamma-ray bursts observed with BATSE that were chosen to investigate the dependence of T, 0 and T 50 on various observed 
profile shapes, (a) Burst 3 B 940217. (fc) Burst IB 910626. (c) Burst IB 910627. (d) Burst IB 91 1 1 18(c) Burst 3 SB 940206. (/) Burst 1B910602. 


The first burst examined, 3 B 9402 17 (Fig. 12a) has a long, T 50 . The measured values of T go and T s0 for the 

complex time profile, with some weak emission near the S/N = 5.5 er case are underestimated; however, we find no 

beginning and end of the event. Figure 14 shows each of the systematic trend in the measured values as a function of the 

reduced time profiles of 3B 940217, except for the case of burst S/N. Note that the size of the statistical uncertainties 

S/N = 200 a. Figures 15a and 15 b show the measured ST go and ST 50 are inversely proportional to the burst S/N. 

values of T go and T 50 as a function of the reduced burst S/N. This was true for each of the 1 2 bursts examined. 

The dashed lines indicate the unreduced values of T go and The profile of the second burst, IB 910626 (Fig. 126), 
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Fig. 13. — Time profiles of six gamma-ray bursts observed with BATSE that were chosen to investigate the dependence of T 90 and T so on various observed 
profile shapes, (a) Burst IB 910503. (fo) Burst 1 B 910601. (c) Burst IB 920517. (tf) Burst 1 B 940429. (e) Burst 3B 9408 17. (/) Burst 1B921207. 


consisted of a short single pulse. We find no systematic the fact that for short duration bursts, with extremely low 

dependence of the measured values of T 90 or T 50 on the values of S/N, the total counts observed from fluctuations in 

burst S/N. In fact, the measured values of T go and T so for background is comparable to the total counts observed 

the reduced bursts were the same as the unreduced value from the burst itself. Thus it is difficult to accurately deter- 

except for the case of S/N = 5.5 <r, where the difference in mine the zero-fluence and total-fluence levels. 

T 90 and r 50 were ~ 140% and ~50%, respectively (both Burst 3B 940206 (Fig. 12e) shows a pulse with substruc- 

still within 1 <r). These large relative differences were due to ture, surrounded by weak extended emission on both sides. 
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TABLE 3 

Results of S/N Reduction of Observed Bursts 


Burst 

Name 


3B 940217 
IB 910626 
IB 910627 
IB 911118 
3B 940206 
IB 910602. 
IB 910503. 
IB 910601 . 
3B 920517. 
3B 940429 . 
3B 940817. 
3B 921207. 


S/N 

(a) 


T><> 

Dr/D a 

(%/<*) 


350,* 200, 100, 50, 30. 15, 5.5 
240,“ 100, 50, 15, 5.5 
84,“ 50, 30, 15, 5.5 
270,“ 200, 100, 50, 15, 5.5 
210,“ 100, 50, 30, 15, 5.5 
21,* 15, 10, 5.5 

590,* 500, 300, 200, 100, 50, 30, 15, 5 5 

406.“ 200. 100, 50, 30, 15, 5.5 

583,“ 400, 200. 100, 50, 30, 15 5 5 

150,“ 100, 50, 30, 15, 5.5 

106,“ 50, 30, 15, 5.5 

438,“ 200, 100. 50, 30, 15, 5.5 


<4/3 

<1(140)71 

<40(70)72(4)* 

<3/1 

<38/ c 

<4/1 

<3/1 

<7/2 

<60/3 

<17/1 

<15(40)73 

< 14(22)72 


* This is the S/N of the unreduced observed gamma-ray burst. 
c The S/N = 5.5 a case had D A > 5 a and D„ x 30%. 

The measured values are more than 5 o from the unreduced value. 
The S/N = 5.5 a case had D R % 140% and D A < 1 a. 

' The S/N = 5.5 a case had D R * 50% and D A < 1 a. 

The S/N = 30 (Tease had D R as 70% and D t < 2a. 

* The S/N = 30 a case had D A < 4 a and D„ as 40%. 

The S/N = 10 o case had D„ * 9% and D A < 1 a. 

‘ The S/N =15 a case had D„ « 280% and D A < 3 o\ 

‘ The S/N = 5.5 a case had D„ s- 10% and D t < 1 cr. 

The S/N = 5.5 a case had D K * 40% and D t <3 a. 

The S/N = 5.5 a case had D„ as 22% and D A < 1 a. 


T,o 

d„/d a 

(%/(t) 


< 3(30) b /2 c 

< 1(50)71 
< 10/2 
<3/1 
<25/ c 
<4(9)71 
<3(280)71/(3)’ 
<3(10)72 
<25/2 

<3/1 

< 10/2 

<5/1 


The differences in T go and T so (see Figs. 15c and 15d) are all 
less than 38% and 25%, respectively. All of the measured 
values became systematically smaller as the burst S/N was 
decreased, similar to the behavior exhibited by the simu- 
lated FREDs in § 3.1.2. 

Burst 3B 920517 (Fig. 13c) contains a weak pulse at the 
trigger time, followed ~4 s later by a strong multipeaked 
pulse of ~ 10 s duration, and then a substantially weaker 
pulse ~15 s later. Though the measured values were not 
systematically larger or smaller as a function of burst S/N, 
we did find that the differences in T go were approximately 
45% for both the 100 a and 50 a cases. The measured values 
of T 50 for these two cases, however, were identical to that 
measured for the original unreduced burst. Large deviations 
at such high values of S/N were not common. Excluding 
bursts 3B 920517 and 3B 940206, and the 30 <r case of burst 
1 B 910627, the relative differences for bursts with S/N > 15 
a were less than 7% and 5% for T go and T so , respectively 
The remaining bursts showed no systematic dependence 
on the burst S/N, other than an increase in the scatter of the 
measured values about the unreduced values. 


4. T 90 AND T 50 DISTRIBUTIONS 

Though the measured values of T 90 (T }0 ) for an event 
observed simultaneously with two different instruments 
may be consistent, the distributions of these parameters 
observed by each instrument may be quite different. This 
can be attributed to different (though likely overlapping) 

populations of burst sources being sampled. 

A distinction must be made between four possible T go 
distributions. The first, and most fundamental, is the parent 
distributions of the intrinsic T 90 s. This is the distribution 
that would be measured by a set of perfect (free of statistical 
noise and background, and having a Detector Response 
Matrix equal to the Identity matrix) bolometric detectors, 


with one detector located in an inertial frame of each burst 
source in the universe. These T go s directly result from the 
burst emission mechanisms. If burst sources are at distances 
of ~1 Gpc, where cosmological effects become important, 
then this distribution may evolve as a function of cosmic 
look-back time. 

The second distribution is the parent distribution of 
extrinsic T go s, and consists of one extrinsic T go for every 
burst that occurs in the universe. These T go s are extrinsic in 
that they are a convolution of the true distribution of intrin- 
sic T go s with the effects of propagation of the radiation from 
the source to the detector. These propagation effects would 
include time dilation and redshift effects resulting from 
cosmological expansion of the universe and any interaction 
of the burst radiation with intervening matter or the back- 
ground radiation. The parent distribution of extrinsic T go s 
would be measured by a perfect bolometric detector 
located in the current epoch, that detects every burst that 
occurs in the universe. 

The third distribution is the selected distribution of 
extrinsic T go s. It is a convolution of the parent distribution 
of extrinsic T go s with the detection criteria and instrumen- 
tal response of the observing detector. The BATSE obser- 
vation of this distribution is shown in Figure 16a this 
distribution will be different for different instruments. 

The fourth distribution is the selected distribution of 
intrinsic T go s. This distribution is obtained by deconvolving 
the propagation effects (but not the detection selection 
criteria) from the selected distribution of extrinsic T go s. 
This is the most fundamental distribution that can be hoped 
to be obtained from the observational data. To obtain the 
parent distribution of intrinsic T go s would require full 
knowledge of the parent distribution of extrinsic T go s 
which is forbidden by the detection criteria of our imperfect 
detectors. At best, the fourth distribution may place con- 
straints on the parent distribution of intrinsic T go s. 
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Fig 14— Time profiles of burst 3B 940217 with various peak S/N ratios, (a) Observed S/N 
Reduced S/N = 30 a (e) Reduced S/N = 1 5 a. (/) Reduced S/N = 5.5 a. 


= 350 a. (b) Reduced S/N = 100 a. (d Reduced S/N = 50 a. id) 


There are a number of factors that contribute to the 
shape of the selected distribution of extrinsic Tg 0 s, resulting 
in different sampling regions of the parent distribution of 
extrinsic T go s for two different instruments. One important 
factor would be a difference in on board trigger criteria. For 
instance, different trigger timescales sample different dis- 
tance ranges. The magnitude of this effect depends on both 
the duration of the burst and on the shape of the time 


profile during the trigger time interval. The shorter trigger 
timescales sample short bursts out to further distances than 
the longer trigger timescales, and the longer trigger time- 
scales sample longer bursts out to further distances. Thus, if 
possible, T 90 (T 50 ) distributions should be examined inde- 
pendently for each trigger timescale, with each distribution 
containing only those bursts that would have triggered that 
timescale at some time. The distributions for all trigger 











1 10 100 1 10 100 
S/N (a) S/N (or) 

Fig. 15. (a)-(d) The measured values of T 90 and T 50 for two of the profiles which have had their S/N reduced. The dashed lines indicate the values of T 90 
and T 50 as measured for the originally observed (unreduced) bursts, (a) Burst 3B 940217. (r) Burst 3B 940206. 






T 90 (Seconds) 


T 50 (Seconds) 


Fig. 16.— (a)-(b) The observed distributions of T w and T so as measured for the BATSE 2B Catalog, (c)-(d) The observed distributions of T, 0 and T 50 for 
those 267 bursts in the BATSE 2B catalog that would have triggered on the BATSE 64 ms trigger timescale. 
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timescales cannot be combined without detailed knowledge 
of the distribution of profile shapes during the trigger time 
interval 

One should not make the mistake of thinking that bursts 
with durations less than t Ug that have low values of peak 
intensity are missing from the selected distributions. This is 
because the term “ peak intensity ” is meaningless unless the 
timescale on which it is measured is also specified. Rather, 
the selected distributions of T go (T 50 ) are flux-complete on 
each individual trigger timescale. Instrument trigger criteria 
require the accumulation of a minimum number of counts 
(C min ), trg on a specific trigger timescale f lrg in a specified 
energy range. All events that result in at least this minimum 
number of counts observed during the trigger time interval 
wilt trigger the instrument. Thus the T go (T 50 ) distributions 
are complete down to a flux of ( C m J t counts in t XT% s. This 
completeness will, of course, shape tne short end of the 
selected distributions of extrinsic T 90 s (T 50 s). The number 
of counts (C min ) lug is calculated based on an assumed 
number of background counts expected during the trigger 
time interval. If the burst duration is less than t uv then this 
assumed background level will be an overestimate on the 
timescale of the duration of the burst. On this timescale, the 
intensity will need to be higher than the threshold intensity 
as measured on the longer timescale t (rg . Thus the region of 
the parent distribution of extrinsic T go s with durations less 
than f lrg will not be completely sampled by any real detec- 
tor, only the population of these bursts with (C min ) tug counts 
on the trigger timescale t lrg will be included in the selected 
distribution of extrinsic T g0 s. 

We stress that this completeness effect should not be mis- 
interpreted as triggering on fluence. The instrument trigger 
criteria are not a function of burst duration; it always trig- 
gers on the number of counts detected above background 
during a trigger time interval. For those bursts with dura- 
tions less than r trg , all source counts (i.e,, the burst fluence) 
may be observed within the trigger interval. 

The trigger threshold of an instrument will also signifi- 
cantly affect the selected distributions of extrinsic T go s. The 
threshold is partly a function of the local background 
environment. For example, the observed local background 
of a detector flown in a near-Earth orbit is quite different 
from that of an interplanetary probe. The threshold is also 
dependent on the collecting area of the detector as well as 
the onboard software triggering criteria. An instrument 
with a lower threshold will observe a sample of less lumi- 
nous bursts that would not trigger a detector of higher 
threshold. In addition, the instrument with a higher thresh- 
old will sample a smaller region of the short-end of the 
parent distribution of extrinsic T go s because of the flux- 
completeness of the selected sample. 

Another trigger criterion often different between various 
detectors is the trigger energy range. Because of the diver- 
sity of the photon spectra of GRBs, triggers examining dif- 
ferent regions of the spectrum will select different samples of 
events (of course with some overlap). 

4.1. Measurements Made with BATSE 

The selected distributions of 433 extrinsic T go s and T 50 s, 
as measured for the set of GRBs in the BATSE 2B catalog 
(Meegan et al. 1994), are shown in Figures 16a and 166. 
Kouveliotou et al. (1993) reported that the distributions are 
bimodal, with a short group (T go < 2 s) and a long group 
(T 90 > 2 s). This result is consistent with previous measure- 


ments of the observed GRB duration distributions (Hurley, 
1992 and references therein). However, the explanation 
behind this bimodality remains unknown. The deficiency of 
bursts is in the region T go % 2 s, which is a factor of ~2 
larger than the largest BATSE trigger timescale of 1.024 s. 
The turnover in the T 90 distribution at values lower than 1 
s, at least in part, reflects the above discussed effect of flux- 
completeness of the selected distributions.The turnover at 
large values of T g0 is most likely real, though the largest 
observed value is probably a result of the finite sampling of 
the parent distribution. 

We have performed a nonlinear least-squares fit to the 
observed T go distribution using a function that represents 
the sum of two log-normal Gaussians, given by 

C(rr v As [ flog T go - log fi s \ 

/<7i °) “^7 p L _ V ^ J 

+ ^ r_/!og7 k -iog 

0L L \ °L 

We find that the best-fit values are n s = 0.60 ± 0.1 1 s, o a = 
1.21 ± 0.14, and /i L = 32.1 ± 2.6 s, a L = 1.06 ± 0.07. The 
reduced y 2 of the fit is 0.816 for 15 degrees of freedom v 
(v = 15), giving P(x 2 , v) = 0.66. 

If one examines only the bursts that would have triggered 
the 64 ms trigger timescale, the bimodality becomes more 
distinct. Figures 16c and 16 d show the T go and T i0 distribu- 
tions for those 267 BATSE bursts that would have triggered 
on at least the 64 ms trigger timescale. For this sample we 
find n, = 0.47 ± 0.07 s, a s = 1.06 ± 0.12, and n L = 34.9 
± 3.2 s, a L = 1.01 ± 0.08. The reduced y 2 of the fit is 0 964 
for 14 degrees of freedom v (v = 14), giving P{y 2 , v) = 0.49. 
Thus the effect of limiting the sample to the smallest trigger 
timescale is to shift the centroid of the short population to a 
lower value of T go and to shift that of the long population to 
a higher value. However, the best-fit values for the single 
timescale sample are within 1 a of the sample containing 
bursts that triggered on any of the three timescales. 

The distributions shown in Figures 16c and 16 d are 
slightly flux-biased. This is because once the BATSE 
onboard burst trigger is set, the 256 and 1024 ms trigger 
timescales are disabled, and the trigger threshold of the 64 
ms trigger timescale is adjusted such that only an event 
stronger than the already triggered burst can retrigger the 
instrument. The three trigger thresholds are reset to their 
nominal values after the readout of burst data from the 
triggered burst is completed. Thus, bursts which occur 
during the readout of data from a previous burst, have a 
peak flux larger than nominal 5.5 a above background on 
the 64 ms timescale, but are weaker than the previous burst 
on this timescale, will not be included in the sample. 

One can examine the region of the selected distributions 
that is nearly free of the effects of flux-completeness, by 
limiting the data to values of T go greater than the shortest 
trigger timescale. Fitting only data with T g0 >0.1 s to the 
model given in equation (35), we find that the best-fit 
parameters are given by // s = 0.56 ± 0.17 s, <x s = 1.68 
±0.57, n L = 34.3 ±3.7 s, and a L = 1.02 ± 0.08. The 
reduced y 2 of the fit is 0.84 for 12 degrees of freedom v 
(v = 12), giving P(y 2 , v) = 0.61. Again, the best-fit values for 
this limited sample are within 1 a of the unlimited sample. 

If GRB sources are at distances of several Gpc, then the 
temporal properties of their emitted radiation should 
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exhibit the effects of time dilation (Weinberg 1972). ( T 90 ) o , 
the value of T 90 for a burst at a cosmological redshift of z as 
observed by a detector at the present epoch, is related to 
(T 90 )„ the intrinsic value of T go that would be observed by a 
detector in a comoving frame near the source, through the 
relation (T 90 )„ = (1 + zX^Wi- Thus, one may predict that 
weak bursts, which are on average further away, should on 
average have larger values of (T 90 ) o . The same effect would 
be expected from an anticorrelation between the peak lumi- 
nosity of a burst and its duration (which would be, for 
example, required if the total energy associated with each 
burst is constant). We have shown that in most cases T 90 
(T 50 ) is not subject to systematic effects that are a function 
of S/N; the exception (primarily due to inadequate back- 
ground modeling) being the case of bursts with weak 
extended emission. These types of profiles show a tendency 
for the values of T 90 (T 50 ) to be systematically underesti- 
mated for low values of S/N. Notice that this works to 
weaken the signatures of cosmological time dilation or of a 
luminosity-duration anticorrelation if either is present in 
the BATSE T go (T so ) distributions, but perhaps more 
importantly, would not introduce either of them. 

4.2. Measurements Made with Other Instruments 

Kargatis et al. (1994) determined T 50 for 171 GRBs 
observed with the SIGNE experiments on board the Venera 
13 and 14 satellites. SIGNE had a smaller detector area 
than a BATSE LAD, required an 8 a fluctuation on a 1.024 
s timescale to trigger, and had a memory size that allowed 
data collection for 62 s. The extrinsic T so distribution 
observed with SIGNE shows a single peak at ~ 6 s, with a 
cutoff at 7^ 0 > 64 s due to the limited memory size, and a 
decaying tail extending from the peak down to small values 
of T so . Below the 64 s cutoff, this peak is consistent with the 
high peak observed in the BATSE T so distribution. The 
SIGNE distribution does not show the second peak at small 
values of T so that is observed in the BATSE distribution 
(Kouveliotou et al. 1993). Kargatis et al. (1994) report that 
the fraction of SIGNE bursts with T so < 0.6 s is 15%, as 
compared to 25% for the BATSE IB catalog. Though these 
numbers are correct, they should have compared the 
SIGNE sample of bursts with those bursts observed with 
BATSE that would have triggered only the 1.024 s trigger 
timescale. In this case, the fraction of BATSE bursts with 
T so < 0.6 s is ~ 15% for both the BATSE IB and 2B cata- 
logs. 

Terekhov et al. (1994, 1995) determined the extrinsic T go 
for 118 GRBs observed with the PHEBUS instrument on 
board the GRAN AT satelite. PHEBUS consisted of six 
bismuth germanate (BGO) scintillation detectors, one on 
each side of the spacecraft. The maximum effective geo- 
metrical area is ~420 cm 2 , ~5 times smaller than that of a 
BATSE LAD. PHEBUS required an 8 o fluctuation above 
background on either a 0.25 or 1.024 s timescale to trigger. 
PHEBUS also triggers in a higher energy range (0.1-1 
MeV) than BATSE (50-300 keV), thus PHEBUS will selec- 
tively choose a sample of bursts from the true distribution 
of observed T go s with harder spectra than those selected 
with BATSE. 

The extrinsix T go distribution observed with PHEBUS 
consisted of a single pulse peaking at ~ 10-20 s, with an 
upper cutoff at T go » 200 s, and a decaying tail extending 
from the peak down to smaller values of T go . The PHEBUS 
distribution does not show a second peak at small values of 


T go . Nevertheless, Terekhov et al. (1994, 1995) report that 
the fraction of PHEBUS bursts with T go < 2 s is 23%, con- 
sistent with that reported by BATSE. This lack of a second 
peak in the PHEBUS sample of 1 18 events may simply be 
due to the better statistics afforded by the larger BATSE 
sample. It may also be due in part to the differing selection 
criteria of the two instruments. Convolving the differing 
trigger energy ranges with the observation that bursts with 
7^o > 2 s show softer spectra in the BATSE data than those 
events with T go < 2 s (Dezalay et al. 1992; Kouveliotou et 
al. 1993) may lead to a deficiency of weak events in the 
PHEBUS data set with T go > 2 s. In addition, the flux- 
completeness caused by the lack of a 64 ms trigger timescale 
on PHEBUS necessitates that PHEBUS samples a different 
set of bursts at the shortend of the extrinsic duration dis- 
tribution than does BATSE. 

5. SUMMARY 

We have provided an extensive description of how T go 
(T J0 ), and their statistical uncertainties, are measured for 
GRBs observed with BATSE. Using simulated simple 
noiseless time profiles, we were able to show that the stan- 
dard software used by the BATSE team to calculate T go 
(T 50 ) is working properly. The only noticeable error, which 
is unavoidable, is that due to the finite binning inherent to 
the data. We also find that in the noise-free cases with 
slowly varying source flux, it can be difficult to distinguish 
source flux from background variations, result in errors in 
T 90 up to «20% (see the case of the FRED profiles with Tj 
of 70s). 

We have investigated the systematic effects associated 
with the measurement of T go and T 50 for simulated bursts, 
to which statistical noise was added. Three profile shapes 
were simulated for a number of combinations of S/N and 
duration. The values of T go ( T so ) measured for the rectangu- 
lar and triangular profiles agreed fairly well with the values 
calculated analytically; we found no evidence for any sys- 
tematic dependence of either T go or T so on the burst S/N. 
Statistical fluctuations in background can have fluences 
comparable to the fluence of short bursts, and become a 
limiting factor in the accuracy of T go (T so ). As in the noise- 
free simulations, the largest errors in the measured values of 
Tgo (Tso) were observed with the case of FRED profiles. We 
found that the measured values of T go ( T so ) were on average 
underestimated by <20% (10%) for the two FRED profiles 
with T d of 30 and 70 s. In addition, we found a systematic 
dependence of T go (T so ) on the burst S/N for the z d = 70 s 
case, and attribute this to inadequate background model- 
ing. For all three profile shapes, the size of the uncertainties 
ST go and ST so were inversely proportional to the burst S/N. 

We examined the measured values of T go (T 50 ) for the 
ensemble of simulated weak (S/N = 5.5 <r) time profiles. A 
systematic dependence of T go (T J0 ) might reveal itself as a 
tendency of either measured value to be larger or shorter 
than the analytic value, The results were consistent with 
some degree of dependence of T go (T 50 ) on burst S/N for 
some profile shapes. 

We studied the effect of including weak source flux in the 
user-defined background intervals. Using our worst known 
case (FRED profiles with r d of 70 s), we found that the 
nearer the background intervals were to the source emis- 
sion, the larger the deviation of the measured value of T go 
(T$o) was from those calculated analytically. The degree to 
which the measurements were affected was essentially iden- 
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tical for both strong (S/N = 50 a) and weak (S/N = 5.5 oj 
profiles. However, it was observed that the background 
intervals could also be too far removed from the event, thus 
not providing enough constraint on the polynomial fit. This 
is particularly true when using higher order polynomials 
to fit small background intervals (limited due to the finite 
time coverage of a data file) separated by large intervals of 
source flux. This will produce large systematic errors for 
low fluence events. 

We randomly selected 10 background models from those 
used in the 2B catalog and combined each model with the 
same simulated burst. Using strong and weak realizations 
of two profile shapes (triangular and FRED), we concluded 
that systematic trends in the measured values due to differ- 
ent background models are not observed. 

Finally, we utilized burst profiles observed with BATSE 
to determine the dependence of T go (T 50 ) on the detailed 
shape of the burst time profile. The signal-to-noise ratios of 
these profiles were reduced to various levels, and we exam- 
ined how the measurements changed as a function of S/N 
for each profile shape. As with the simulated time profiles to 
which noise had been added, we found that the size of the 
uncertainties ST go and <5T 50 are inversely proportional to 
the burst S/N. Only one burst (3B 940206), consisting of a 
single pulse surrounded by weak extended emission on both 
sides, exhibited a systematic behavior; this behavior was 


similar to that observed for the simulated FREDs (see Figs. 
7c and Id) in which the source flux varied on the same 
timescale as the background. The other selected bursts 
showed no evidence for a systematic dependence of T go 
(T 50 ) on the burst S/N. In addition, we found the changes in 
Tgo (T so ) induced by a reduction in S/N depended on the 
detailed shape of the time profile near the times r f ; in some 
cases, T go was less sensitive to a reduction of the S/N than 

*We have distinguished between parent and selected dis- 
tributions of intrinsic and extrinsic T go s (T 5 q s). The only 
distribution that is observed is the selected distribution of 
extrinsic T go s. Some of the factors that shape this distribu- 
tion are discussed. The selected distribution of extrinsic T go 
s as observed with BATSE was presented and compared to 
previously reported observations made with other GRB 
experiments. The deficiency of bursts at T g0 ~ 2 s in the 
BATSE-selected distribution of extrinsic T go s is well above 
the instrument’s largest trigger timescale (1.024 s) and it is 
unlikely a result of the effects of flux-completeness. 
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APPENDIX A 


BATSE INSTRUMENTATION AND DATA TYPES 

BATSE consists of eight detector modules mounted on the eight corners of the CGRO spacecraft in such a way as to 
provide nea/aH-sky 0 coverage. There are two gamma-ra, de, errors within each i mod ale, a Large Area Detector (LAD) and a 
Spectroscopy Detector (SD). The LAD, an uncollimated Nal(Tl) scintillation detector with a surface area of ~2025 cm , 
optimized for sensitivity over the energy range -20 keV-2 MeV with modest energy resolution. This study does not use SD 

ILVrSE so r Is th e ^o bs e r v ed countsTnTo fou^discriminator energy channels every 64 ms. The energy ranges of these four 
channels are approximately 20-50 keV, 50-100 keV, 100-300 keV, and the last channel ,s an , SXSr 
threshold of -300 keV. The discriminator data, summed over the energy range 50-300 keV (channels 2 and 3) 
continuously monitored for the satisfaction of a set of criteria, upon which a burst trigger is dec ared These criteria 'nclude a 
5 5 a increase above background simultaneously in two or more detectors, on at least one of three timescales (64 256 and 
1024 ms). When these criteria are satisfied, the instrument enters the burst acquisition mode whereby various burst data types 

ar Th?d™iU^ (PREB), the discriminator science (DISCSC), the large area discriminator 

(D ]JcL“dTa"nd1he n ,, m e-,»gid even.V, data. Each , of These da, a (ypcs fTSfi 

energy channels. The PREB data cover the 2.048 s prior to the burst trigger time T lrg . The DI £CSC data st^ 
for - 240 s prior to 1992 December 17, and for - 550 s thereafter. The time resolution of the PREB and DISCSC data types is 
64 ms. The DISCLA data are always available, except during times of CGRO telemetry gaps, and have a time resolutio 

1 While not in a burst acquisition or data readout mode, the TTE data (consisting of the time of arrival and energy channel >f 
each interacting photon) are continuously being recorded in a ring buffer memory of fixed size. Once a burst trigger >sset,the 
most recent i of the TTE memory are stored and the remaining memory is filled as the data are accumulated. Thus, the sta 
time of the TTE data is dependent on the burst intensity just prior to the burst trigger time, andthe end me of the TTE 
depends on how fast the memory is filled after the trigger time. The intrinsic time resolution of TTE data is 2 (is, but these data 

ar AsTi^the 0 DISclTdata °due f to^elemefry ^apT the'pR EB, DISCSC, and TTE data types may not be available for all 
bursts In addition TTE data are not available for those bursts that are overwrites. An overwrite burst is a burst that triggers 
the instrument while it is still in the burst acquisition or data readout mode. Once the burst trigger is set, the 64 ™ t r igge 
hreshotrs adjusted such that only a burst stronger than the current trigger is allowed to retrigger the instrument. The 256 
and 1024 ms thresholds are disabled during burst acquisition and the subsequent data readout. Once the burst trigger isset 
■ t . TTF mpmnrv are filled the writing of TTE data to memory is halted and not restarted until after the data readout is 
over Therefore the TTE data associated with an overwrite burst (occurring before the accumulation of TTE is restarted) wi 
belong to the previous burst. Lack of TTE data affects only short duration bursts whose time histones are not well resolved 
with the 64 ms data resulting in large relative errors in the T go ( T $0 ) measurements. 
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ANALYTIC CALCULATIONS OF T 90 AND T 50 FOR THREE SIMPLE TIME PROFILES 


The simple forms of equations (20), (23), and (24) allow the values of 7^ 0 and T 50 to be calculated analytically using 
equations (2), (3), and (4). For rectangular bursts, represented by equation (23), equation (2) becomes 


/ _ &Pdt 

100 Pdt ’ 


(36) 


where P is the peak source count rate summed over all energy channels. Solving for x f gives 


T / = 



(37) 


where the duration T d — x e — z s . Using equations (3) and (4) we find that for simple rectangular burst profiles, 

T 9O = 0.9T d and T 50 = 0.5T, . (38) 

The analytic values of T go (T so ) for the FRED bursts, which are represented by equation (24), can also be calculated 
straightforwardly. Equation (2) is solved for z f , giving 

x f = x s -x d ln(l -//100). (39) 


Equations (3) and (4) result in 


7,0 = ^ In (95/5) and T 50 = z d In (75/25) . 


(40) 


Analytically calculating the values for T go (T 50 ) of triangular bursts, represented by equation (20), is slightly complicated by 
the discontinuity at the peak time i pk . Using equation (1), the total source counts S tol is given by 

S.ot = l( m - ") T pk - l m ?s + i nz e +(b- d)z pk + z e d - z s b , (41) 

where m, n, b, and d are given by equations (21) and (22). For a simple triangular burst, the exact analytic form ofz f depends 
on the size off relative to the fraction / pk of the total source counts S tot that have been observed by the time t = r pk . 

Iff < f then z f < r pk and equation (2) becomes 

J_ = ]\{ (mt + b)dt 
100 S lot 


This gives a quadratic equation for x f , 


where g is a constant given by 


jmzj + bzf-g = 0 , 


9=\ mzf + bz s + ^ S lot 


(43) 

(44) 


The correct root t f of the quadratic equation is determined by a set of physical constraints on z s . This equation is valid for 
/</ k , giving the first constraints z s <z f < r pk . Additional constraints are provided by the set of conditions z 5 < r 25 < 

T 75 T 95* 

If/ > f pk then T f > r pk and equation (2) becomes 

/ = jy {mt + b)dt + (nt + d)dt 
100 S lot 


This gives a different quadratic equation for t f , 

jmj + z f d + h = 0, (46) 

where h is a constant given by 

h = jm(z 2 k - zf) + b(z pk - zf - ±nzf k - z pk d - S 101 . (47) 

The correct root t f is again determined by a set of physical constraints on x f . This equation is valid for / >/ pk , giving the first 
constraint t pk < x f < x e . The additional constraints are again provided by the set of conditions t 5 < r 25 < t 75 < t 95 . 

T go and T so are calculated as before using equations (3) and (4), and the equations (43) and (46) for x f . Because of the 
dependence of the form of x f on its position relative to r pk , no simplified expression of either T 90 (T 50 ) can be obtained for the 
simple triangular bursts. 
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